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Abstract 

We discuss the potential of high sensitivity mid-IR and far-IR spectroscopy to proof the 
physical properties of active nuclei and starburst regions of local and distant galaxies. For local 
galaxies, it will be possible to map the discs and ISM through the low ionization ionic lines and 
a variety of molecular tracers, such as OH, H2O and high-J CO. At increasing distance, most of 
the ionic nebular lines (typical of stars and AGNs) are shifted into the FIR, making possible to 
compare the observed spectra with those predicted by different evolutionary scenarios. At the 
very high redshift of 10-15, sensitive mid-to-far-IR spectrometers, such as those planned to be 
at the focal plane of the future SPICA mision, could be adequate to detect the H recombination 
lines excited in the HII regions around population III stars, if these stars happened to reside 
in large clusters of more than 10 5 members. 

1 Introduction 

The spectroscopic observations of the Infrared Space Observatory (ISO) opened a new window 
for the study of the physical and chemical properties of IR-bright, ultraluminous infrared galaxies 
(ULIRGs) and active galactic nuclei (AGN). Most of our present knowledge of the mid-to-far 
infrared spectra of active and starburst galaxies is derived from ISO spectroscopy (Fischer et al 
1999, Sturm et al. 2002, Verma et al. 2003). The mid-IR spectral range includes most of the 
fine-structure lines excited by the hard radiation produced by black hole accretion as well as those 
mainly excited by stellar ionization (e.g. Spinoglio & Malkan 1992) and thus represents an essential 
tool to distinguish between the two processes, especially in obscured nuclei suffering severe dust 
extinction. In a complementary way, the far-IR range contains a large variety of molecular (OH, 
H2O, high-J CO) and low excitation ionic/atomic transitions, in emission or in absorption, that 
can reveal the geometry and morphology of the circumnuclear and nuclear regions in galaxies. In 
particular these latter could trace the expected conditions of X-UV illuminated dusty tori (see e.g. 
Krolik & Lepp 1989 for predictions of high-J CO lines), whose presence in type 2 active galaxies 
is foreseen to reconcile the typel/type2 dichotomy. 

The far-IR spectra of local IR-bright and ULIRGs galaxies, as measured by ISO-LWS (Fischer 
et al.1999), showed an unexpected sequence of features, as can be seen in Fig.l, from strong 
[OIII]52, 88/im and [NIII]57/^m line emission to detection of only faint [CII]157/im line emission 
and [OI]63/im in absorption. The [CII]157/im line in 15 ULIRGs (Lju > 10 12 Lq) revealed an order 
of magnitude deficit compared to normal and starburst galaxies relative to the FIR continuum. 
Non-PDR components, such as dust-bounded photoionization regions, generating much of the FIR 
continuum but not contributing significant [CII] emission, can explain the [CII] deficiency. Such 
environments may also explain the suppression of FIR fine-structure emission from ionized gas and 
PAHs, and the warmer FIR colors of ULIRGs (Luhman et al 2003). 

Examples of far-infrared molecular lines detected by ISO include: OH at 34.6/im, 53.3/im, 
79.1/zm and 119/xm connecting to the ground state, giving rise to absorption in most cases; OH 
at 98.7^m, 163/mi at higher energies, in emission; H2O at 73.5/im, 90^m, 101/im, 107/im, 180/im 
in emission or absorption. From the absorption lines the gas column and the abundance can be 
derived, from the emission lines we can measure the gas cooling and constrain the temperature 
and density of the warm (50<T<500K) molecular gas. 



LWS observations of Arp 220 show absorption in molecular lines of OH, H2O, CH, NH, and 
NH3, as well as in the [0I]63/J,m line and faint emission in the [CII]158/im line. The molecular 
absorption in the nuclear region is characterized by high excitation due to high infrared radiation 
density (Gonzalez- Alfonso et al 2004). The LWS spectrum of the prototype Seyfert 2 galaxy 
NGC1068, beside the expected ionic fine structure emission lines, shows the 79, 119 and 163/im 
OH rotational lines in emission, not in absorption as in every other galaxy yet observed. Modeling 
the OH lines suggests the emitting gas lies in small (O.lpc) and dense clouds (~ 10 4 cm -3 ) most 
probably in the nuclear region and is potentially a signature of the torus (Spinoglio et al 2005). 

From the pioneering work of ISO it appears that both the mid-infrared and the far-infrared 
are powerful tools to study the physical conditions in galaxies, often obscured by dust. The 
future far-infrared spectrometer ESI (Swinyard 2006) at the focal plane of the helium cooled 
3.5m class telescope of the SPICA mission (Nakagawa 2004, Onaka & Nakagawa 2005), due to 
its high sensitivity, will be able to detect FIR lines in local galaxies, near-to-mid infrared lines 
and dust/PAH features in distant galaxies and possibly reveal the presence of large clusters of 
population III stars through H recombination lines. 
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Figure 1: The ISO-LWS spectra of six IR- bright galaxies (adapted from Fischer et al 1999). The spectra are 
shifted and ordered vertically according to apparent excitation (not in order of relative luminosity or brightness). 
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2 Infrared line diagnostics through fine- structure emission 
lines 



Mid-IR and far-IR spectroscopy of fine-structure emission lines are powerful tools to understand 
the physical conditions in galaxies from the local universe to distant cosmological objects. Figure 
2 shows the critical density (i.e. the density for which the rates of collisional and radiative de- 
excitation are equal) of each line as a function of the ionization potential of its ionic species, in 
three different redshift ranges, one for each frame, for which the rest-frame wavelength of the 
line is shifted in the far-infrared range. This diagram shows how these lines can measure two 
fundamental physical quantities (density and ionization) of the gas. In the figure are shown with 
different symbols the lines originated in different astrophysical emission regions at work in galaxies: 
the photodissociation regions (PDR), the stellar ionization in HH/starburst regions, the AGN 
ionization and the high ionization of coronal type emission. Choosing two lines along the abscissa 
and two along the ordinate, their ratios give, respectively, a good estimate of the ionization and 
the density of the gas in the region (see, e.g., Spinoglio & Malkan 1992). 
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Figure 2: Fine-structure lines occurring in the 30-200/mi range at three different redshift ranges (the rest wave- 
length is reported near each line label). The lines are plotted as a function of the ionization potential and their 
critical density. Different symbols are used for the lines from the photodissociation regions (PDR), shown as squares, 
the stellar/HII region lines, shown as triangles, the AGN lines, depicted with circles and the coronal lines, indicated 
with stars. The two [OI] lines have been plotted - for graphical reasons - at a ionization potential higher than their 
zero value. 

It appears from the figure that in the far-IR range, spectroscopy of the fine structure lines can 
probe in the Local Universe the two regimes of photodissociation (PDR) and stellar ionization. 
The stellar emission (e.g. in starburst galaxies) can be probed up to high z, using many lines in 
the rest-frame spectral range of 3 < X(fim) < 30. The ionization from AGN can be probed in the 
redshift range l<z<5 and the coronal emission regions are probed by near-IR lines shifted into the 
far-IR at a redshift of z ~ 5. 
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3 Fine-structure line detection in Seyfert galaxies at in- 
creasing redshift 

The ISO spectrometers detected in the bright local Seyfert type 2 galaxy, NGC1068, many of the 
lines plotted in Fig. 2 at flux levels of 5-200 x 10~ 16 W m~ 2 (Alexander et al 2000, Spinoglio et al 
2005). This object contains both an active nucleus and a starburst. Considering this galaxy as a 
template object, we computed the line intensities expected at redshifts ranging from 0.1 to 5. For 
simplicity, we adopted an Einstein-De Sitter model Universe, with tt\ = fl vac — and Qm= 1, 
H =75 km s _1 Mpc -1 . The luminosity distances have been derived using: 

d L (z) =2c/H [l + z- y/(l + z)] (1) 

The results are reported in Table 1, where the line intensities are given in W m~ 2 . We have 
assumed that the line luminosities scale as the bolometric luminosity and we have chosen two cases: 

A) a luminosity evolution proportional to the (z+1) 2 , consistent with the Spitzer results at least 
up to redshift z=2 (Perez-Gonzalez et al 2005) (A values in the table); 

B) no luminosity evolution (B values in the table). 

Because the star formation process in galaxies was much more enhanced at z=l-2 than today, 
we consider reasonable to adopt the model with strong evolution (model A) for the stellar/HII 
region lines and the "no evolution" one (model B) for the AGN lines. 

We note that the dependence on different cosmological models is not very strong. The popular 
model with flu— 0.27, fl vac —0.73, H =71 km s _1 Mpc -1 shows greater dilutions, increasing with 
z, by factors of 1.5 for z=0.5 to 2.5 for z=5. 

4 Detection of the first stars 

Primordial stars are born in a Universe which was very poor in metals, for this reason they should 
form in very massive clouds which are expected to originate in very massive clusters populated by a 
large number of massive stars. Following the suggestion by Panagia (2005) to detect the primordial 
stars though their HII regions nebular emission, we modeled the emission of primordial clusters 
using the model stellar atmosphere by Schaerer & de Kotcr (1997) and a standard photoionization 
model. 

We have used a stellar model with mass M=120 M© (the largest value of the mass among 
the models by Schaerer & de Kotcr), effective temperature T e ff— 53,400K, luminosity L=2 x 
10 6 L Q , and radius R=15.6 R©. We computed the spectrum of an HII region originated from 
this star, using the CLOUDY v96 photoionization code with the following parameters: internal 
radius R(int)=0.1pc, external radius R(ext)=2.5pc, metallicity of 10~ 6 Z©, Hydrogen density of 
10 3 cur 3 . 

We found that to obtain flux densities comparable to the spectroscopic detection limit of JWST, 
about 10 5 stars are needed in the cluster. We present in Fig. 3 the synthetic spectrum of such a 
cluster at a redshift of z=15. It can be seen from the figure that such clusters, if they existed in 
the primordial Universe, can be detected through Hydrogen and Helium recombination lines by 
JWST in the mid-IR (e.g. in the Lya, Ha, He I 1.08/xm lines) and by SPICA in the far-IR (e.g. 
in the Paa and Bra lines). The lines that can be observed in the far-IR are better than the Lya 
line, because this latter can be severely absorbed by the dust, which could be present in the HII 
region, if one allows moderate metallicities of the order of 10~ 3 Z© (Panagia 2005). Moreover, the 
Lya line can also be absorbed by HI in the intergalactic medium (Miralda-Escude & Rees 1998). 

The sensitivity of a mid-to-far-IR spectrometer onboard of SPICA should be better than that 
one of JWST, at wavelenths larger than 10 fim, because it will not be limited by the thermal 
background of the mirror. We conclude that the estimated number of members of the cluster can 
well decrease by a factor 10, to the value of 10 4 stars, if we consider stars with a larger mass, e.g. 
of the order of 150-200 M©, for which reliable models are not available. 
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Figure 3: Redshifted (z=15) starburst spectrum of 10 5 stars ionizing HII regions Z = 10 6 Zo. Flux is scaled at 
Earth. The spectral ranges of JWST, ESI onboard of SPICA and ALMA are indicated. The spectroscopic sensitivity 
of JWST is given for t=4xl0 5 sec, S/N=10, R=1000. 
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Table 1: Line intensity predictions in W/m 2 for NGC1068 at different redshifts 
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Notes: 


for each transition 


in each column, the first line 


gives the observed intensity of NGC1068, the other lines 



list for the given redshift the predicted intensity value for models A (top value) and B (bottom value). The blank 
values correspond to observed wavelengths outside the ESI-SPICA range ~30-200/im. 
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